Hydrodynamic Limit for the Spin Dynamics of the Heisenberg Chain 
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We show that Quantum-Monte-Carlo calculations of the dynamic structure factor of the isotropic 
spin-1/2 antiferromagnetic chain at intermediate temperatures corroborate a picture of diffusive spin 
dynamics at finite frequencies in the low-energy long wave-length limit and are in good agreement 
with recent predictions for this by J. Sirker, R. G. Pereira, and I. Afileck ^arXiv:0906.1978t ^l]. 
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PACS numbers: 

The one-dimensional (ID) Heisenberg XX Z antiferro- 
magnet 



H 



where J > is the exchange coupling, S'f'^ are spin-1/2 
operators on site Z, and A is the exchange anisotropy is 
one of the best studied strongly correlated many-body 
system. Its magnetic transport properties however, re- 
main an open issue Spin transport in the Heisenberg 
chain is directly related to carrier transport in ID corre- 
lated spinless fermion systems, via the Jordan- Wigner 
transformation, and therefore is of great interest in a 
broader context. Linear response theory 01 shows the 
zero momentum, frequency dependent spin conductivity 
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to consist of the Drude weight 
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and a regular spectrum 
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where j = jq=o is the z-component of the spin current 
with jg = (iAJ/2) cxp(— i(7Z)(S';~S';'*^^ — S^Sf^^^) and 
m, n are the eigenstates with energies Em,n- 

The Drude weight has been under intense scrutiny for 
more than two decades. However, no generally accepted 
picture has emerged. A nonzero Drude weight would im- 
ply dissipationless transport in a correlated system 
despite the fact that [j, H] ^ for the XX Z model. Here 
we give a brief summary regarding the status of this is- 
sue and refer to [H and refs. therein for a more extensive 
summary. At T = and in the massless regime |A| < 1 
of the XX Z chain, the zero temperature Drude weight 
is known to be finite ji]. At T 7^ 0, Bethe-Ansatz (BA) 
calculations arrive at contradictory results regarding the 



temperature dependence of D(T) 0, 0, 0|- The same 
holds for the question whether D{T > 0) is finite or not 
at the SU{2) symmetric point A = 1 (1, 0|. Recent nu- 
merical studies using QMC [1,0], exact diagonalization 
(ED) at zero 0, [foTflll. [l^ . as well as finite magnetic 
fields [13II, and master equations (li . [l5| are consistent 
with D 7^ for |A| <1 and T > 0, supporting a baUistic 
contribution to the conductivity at finite temperatures. 
Recent time-dependent density-matrix renormalization 
group (tDMRG) studies have given evidence for ballis- 
tic spin dynamics for |A| <1 in the out-of-equilibrium 
case (l6l |. 

The regular finite-frequency contribution cr^egl"^) has 
been considered by ED studies 17, which however 
leave many open issues. Very recently, spin diffusion 
has been conjectured to govern the low-frequency spec- 
trum of the regular conductivity based on real- 
time transfer matrix renormalization group (tTMRG) 
and a perturbative analysis using bosonization. The lat- 
ter provides for an approximate expression for the Fourier 
transform of the retarded spin-susceptibility Xret i<l,t) = 
iQ{t){[S^{t),Slg]), which reads 
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where at A = 1, K ~ \, v ~ 7r/2, 2^b — ''^9^T, 
b = 52/4 _ ^3(3 _ 87rV3)/32 + v^T^/tt, and c = 
— 3g^/32 — a/ST^/tt have been obtained by per- 
turbative expansions (PE) at T ^ J in pow- 
ers of the running coupHng constant l/g + ln{g)/2 = 

ln(yV2exp(G-h l/4)/r) and G w 0.577216 .. . is Eu- 

ler's constant [iot . 

Some remarks are in order. First, for a; ^ 7, eqn. 
(O displays a diffusion pole with a diffusion constant 
r = (1 -I- c)v'^ / {■Kg'^T). I.e. within this approxima- 
tion the spin dynamics of the Heisenberg chain would 
allow for a plain hydrodynamic limit. Second, eqns. 
(HI and ([6]) do not incorporate the finite width of the 
spectral function x" {q-i^) ~ I'^lXretiq,^)]/''^ at T = 
0, which is dominantly set by the two-spinon contin- 



2 



uum. However, at g/vr <^ 1 the latter width is of or- 
der nJq'^/lQ, which for those wave vectors and tem- 
peratures which we will be interested in is negligible 
against 75. Third, for any finite momentum q ^ 0, 
the isothermal susceptibility Xq = d,uJX"iQj^)/^ ob- 
tained from eqn. ([5]) is identical to the isolated sus- 
ceptibility Xiet('?, 0) = J^^dujx"{q,i^)/{uj - iO^), since 
X"{q 7^ 0, w ^ 0) oc w. Therefore Xq = K/{2-kv{1 + c)). 
Furthermore, the isothermal susceptibility of the Heisen- 
berg model is a continuous function of q. Its limiting 
value limq_,o Xq = Xo at zero momentum is known from 
thermodynamic Bethe Ansatz (TBA) [13, HI. Therefore 
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should be satisfied, where xpe is a known PE of the TBA 
result [13, [2l|. This implies, that apart from eqns. (A2) 
and (B3) of ref. also the ratio K/v requires to be 

renormalized off from 2/tt . 

The spectral function x"('Z, w) is related to a'^^g{uj) by 
means of the lattice version of the continuity equation 
dtSq = (1 - exp(-ig)) j, through 



'^'regi^) = 1™ — x" ('7, • 
q^Q q^ 
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Therefore, the spectrum of the regular part of the optical 
conductivity can be deduced from eqns. ^ and ([6]). 

The main goal of this work is to analyze, to which 
extent eqns. ^ and |[6]) are consistent with QMC calcu- 
lations. The significance of such comparison is with the 
regular part of the spin conductivity. It will not clarify 
the size of the Drude weight, as any discrepancy aris- 
ing may be due to partial spectral weight transfer into 
a Drude weight. Furthermore, we focus on the isotropic 
point A = 1, which may be different from the anisotropic 
case. To begin, we note, that eqns. ^ and ([6]) approxi- 
mate the on-shell part of the spectrum for |w ± vq\ ^ T. 
Yet, similar to the comparison with tTMRG in eqns. 
(C2) and (C3) of ref. we will assume them to be valid 
for all uj. Furthermore, Xq is known to monotonously 
increase for the Heisenberg model as q ^ n/2. How- 
ever, Xq = K/{27Tv{l + c)) from bosonization is momen- 
tum independent. Therefore, a momentum dependence 
K Kq, V Vq- albeit weak at g ^ 1 - is to be allowed 
for, when matching up eqns. ([5]) and ^ with QMC. 

We perform the comparison to QMC by transforming 
Xret{q,uj) onto the imaginary time axis 



2 ^ cos(tj„T)x((7, ujn) - x{q^ 0) 
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The main point is, that a corresponding xqmc i<l, t) can 
be obtained directly from QMC, following preceding work 
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Figure 1: Imaginary-time susceptibility XQmc{<1,t) at g = 
7r/64 on 128 sites, for two temperatures T, fitted to x(9i ''') 
from eqn. ^ (lines) in three ways, namely: "Jq^QMC opti- 
mized (solid), 7s taken from ref. |l9l | (dashed), and 7 forced 
to zero (dotted). The index 'i' on the y-axis refers to xil-,''') 
from eqn. ((QJ for the lines in panels a) and b) as well as to 
QMC for the symbols in panel a). Panel a) Global behav- 
ior ot XQMciq,T)/xQMc{q,0) for T/J = 0.1 (QMC, squares) 
and 0.25 (QMC, circles). In this panel the three fits (lines) 
are indistinguishable on the scale of the plot. Panel b) Error 
2a of XQMc{<},f) for each r evaluated (error bars) and dif- 
ference x(9, t) — XQMc{q, t) between QMC and the three fits 
(lines). 2(T for the QMC data is O(10~^). Plots corresponding 
to T/J = 0.1 have been shifted by 2 x 10~^ 



employing the stochastic series expansion method |23l |. 
This involves only the statistical error, which is well con- 
trolled. Uncontrolled sources of error, due to e.g. trans- 
formations to real or Matsubara frequencies, do not oc- 
cur. X (9, t) is gauged against xqmc {q, t) by fitting Kq, 
Vq, and 7g at small momentum, while retaining b and c 
as given by bosonization. This is justified, because the 
latter two constants do not enlarge the space of fitting- 
parameter, as any modification of them can be absorbed 
into a renormalization of Kq, Vq, and 7,. Regarding the 
temperature range, we confine ourselves to T/J < 0.25. 
This is motivated by the PE to 0{g^,T'^) for thermody- 
namic properties to agree rather well with QMC results 
up to T/J < 0.1 jiil, while for T > 0.25 the PE starts 
to fail significantly. 

Fig. [Ushows the result of the comparison of QMC with 
eqn. ^ for the smallest non-zero wave- vector g = 7r/64 
of a 128-site system for two temperatures T/J = 0.1 and 
0.25 allowing for three different choices of 7^, namely (i) 
lq,QMC as optimized by fitting, (ii) 75 taken from the 
bosonization, and finally (iii) 7^ = forced to be zero 
[iil . The upper panel b) of this figure clearly demon- 
strates, that QMC is inconsistent with 7^ = and that 
increasing 7g above zero improves the quality of the fit. 
In particular the best fit, i.e. for 7g,QA/c, is identical 
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Figure 2: Momentum dependence of the renormalized Lut- 
tinger parameter Kq, spinon velocity Vq, and scattering rate 
7i},QA/c/7s for the first non-zero six momenta on a 128 site 
system for two temperatures T/J = 0.1 (white symbols) and 
0.25 (black symbols). Note that 'yq.QMC /ib for T = 0.1 has 
been scaled by 2 to fit into the plot. 



within the standard deviation 2a (error bar) to QMC for 
almost all r g [0, /3] at both temperatures. Yet, we find 
Iq.QMC > 7B, and moreover there are systematic oscilla- 
tory deviations. While the latter seem a subdominant ef- 
fect, which could be due to the on-shell approximation in 
eqns. |[5| and these deviations may also indicate rel- 
evant corrections to diffusion and should be investigated 
in future studies. We emphasize the vertical scale on 
panel b) of fig. [1] which demonstrates that high-precision 
QMC is mandatory for the present analysis. Fig. [1] is a 
central result of this work. It shows that QMC is consis- 
tent with a dynamic structure factor of the isotropic an- 
tiferromagnetic Heisenberg chain which is approximately 
diffusive at intermediate temperatures in the long wave- 
length limit with a diffusion kernel (1 -I- c)w^/ {'ijq^QMc)- 
Any momentum dependence oi ^q,QMc , to be discussed 
later, imphes corrections to this diffusion. Next, and to 
further support our approach, we will also discuss the 
Luttinger parameters we find. 

In table U we compare the parameters obtained from 
the fit to QMC with results from TBA, PE and tTMRG. 
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0.0190 





Table I: Columns 2 and 3: Comparison of Xq.QMC ~ 
KJ{2nvq{l + c)) from QMC at g = 7r/64 with xo from TBA 
|2d | and xpe from the l.h.s. of eqn. ([7]). Columns 4, 5, and 6 
display 7 from bosonization, tTMRG, and QMC. 



This table shows, that Xq,QMC = Kq/ {2T:Vq{\ -\- c)) at 
q — 7r/64 is in excellent agreement with the isothermal 
susceptibility at g = from the TBA for both tempera- 
tures which we have studied. This result should not be 
confused with the well know agreement between static 
QMC and TBA for the isothermal susceptibility 
but rather it is a satisfying consistency check for our ap- 
proach. In fact, fitting the imaginary-time transform of 
an approximate i-e. eqn. |[5|), to QMC could 

require values for Kq, Vq, and 7g which deviate from ex- 
actly known values for these quantities on a scale which 
is unrelated to the error 2a of the QMC. As will be shown 
later the variation of Kq and Vq with momentum is very 
weak as g ^ 1, i.e. we expect no relevant change for 
Xq,QMC as (7 ^ 0. Yet we are tempted to point out, that 
Xq=Tr/64.QMC in table Uis barely larger than xo, which is 
consistent with the momentum dependence for the exact 
Xq- The fact that Xq,QMC /xpe > 1 and is increasing 
as T increases, evidences that xpe on the l.h.s. of eqn. 
([7]) increasingly underestimates the TBA result as T in- 
creases beyond T/J > 0.1. In fig. [T]we have shown, that 
lq,QMC 7^ IB- Yet, table U demonstrates that jq^QMc 
and 7b are comparable to within factors of order 2. Most 
important, the relaxation rate "fq^QMC we find is much 
larger than the width of the two-spinon continuum, yet, 
very small compared to temperature 7<j,qa/c ^ T. We 
note, that fits to tTMRG O at T/J = 0.2, lead to 
Iitmrg/ib ~ 0.64. 

Next we discuss the momentum dependence. Fig. [2] 
displays all three fit parameters Kq, Vq and jq^QMC ver- 
sus the first six non-zero momenta and the two temper- 
atures T/J = 0.1 and 0.25 which have also been con- 
sidered in fig. [21 Vq and jq^QMC have been normalized 
to their values given by bosonization, i.e. tt/2 and 75. 
Obviously all momentum variations are very smooth and 
rather weak. As can be seen from this figure, most of the 
renormalization of the ratio Kq/vq from its bare value of 
2/tt stems from Kq > 1. The spinon velocity Vq deviates 
slightly from tt/2, however only to within 0(1%). As 
discussed in the previous paragraph, this is necessary to 
obtain an optimum fit of the QMC to the approximation 
eqn. ^ and does noi imply that QMC is at variance with 
the bare spinon velocity. Kq displays a very weak upward 
curvature, while Vq shows a small downward curvature. 
The latter can be understood in terms of the 0{q^) cor- 
rections to the linear on-shell dispersion io{q) which are 
not contained in bosonization. The combined momentum 
dependence oiKq/vq leads to the expected increase of the 
static susceptibility with q. Finally, ^q,QMC /iB also dis- 
plays a weak momentum dependence which is larger for 
T/ J = 0.1. The latter may signal the onset of finite size 
effects. In fact, jq^QMC 7^ implies a length scale I of 
order 0{v / {2'yq^QMc)) for the regular current relaxation. 
I is less than the system size for both temperatures stud- 
ied. Yet, 128/Z w 9 for T/J = 0.25 and and 128// w 3 
for T/J — 0.1. With momentum dependence, 7 as ex- 
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Figure 3: t<-'nO"Qj\/c(9, i^n) from QMC for the first three non- 
zero Matsubara frequencies aj„ = 2nnT and wave vectors q = 
mr/M, with n = 1,2, and 3 as compared to LOnXil^^n) / 
using eqn. ([9]) with 7 = (dashed) and 7 = 'yq,QMC (solid) 
on a 128 site system for a) T / J = 0.1 and b) 0.25. (See text 
regarding statistical error.) 



tracted from a real-space quantity [l9|| will differ from 
that obtained by QMC at fixed small momenta. 

While the preceding has been exact up to the sta- 
tistical error of the QMC, we would like to conclude 
this work by speculating on the line-shape of the reg- 
ular part of the conductivity on the imaginary fre- 
quency axis at w„ = 2TrnT. In principle this requires 
a careful analysis of the error introduced by the Fourier 
transform XQMc{q,^n) = J^^^^ exp{iujnT)xQMc{q,T)dT. 
This error will increase as w„ increases. Here we re- 
frain from analyzing this, since our goal is merely to 
demonstrate to which extend our QMC data discrim- 
inates between a conductivity with 7 = and one 
with 7 — jq.QMC 7^ 0. To this end fig. [3] displays 
uJna-QMciq,^n) = ^nXQMc{q,t^n) / as Compared to 
a;„cr(g,w„) = Lo'^^xiqi^n) / q^ with xiq^^n) taken from 
eqn. ^ and with 7 = or 7 = "fq.QMC- Without any 
further ado, this figure clearly demonstrates that 7 = 
in a{q, tOn) from eqn. ^ and ^ is inconsistent with our 
QMC which however agrees very well with a{q,ujn) for 
7 = lq,QMC |23- This implies that QMC is consistent 
with a Drude type of behavior of the frequency depen- 
dence of the regular conductivity with a relaxation rate 
'^lq,QMC- While future studies, may focus on finite size 
scaling, to perform the limit of g — > 0, as required in eqn. 
([8]), this is beyond the scope of the present analysis. 

In conclusion QMC is consistent with spin dynamics 
of the isotropic ID Heisenberg antiferromagnet which is 
primarily diffusive in the long wave-length limit and at in- 
termediate temperatures, implying a regular part of the 
spin conductivity with a finite relaxation rate 7 ^ T. 



This corroborates recent findings by bosonization and 
tTMRG. Our analysis does not allow conclusions on the 
pending open questions on the Drude weight at A = 1, 
yet based on the numerical evidence for D{T > 0) > 0, 
our findings may open up the intriguing possibility of 
a finite temperature dynamical spin conductivity of the 
isotropic Heisenberg model which comprises of both, a 
finite Drude weight and a regular part with a very large 
mean free path at low temperatures. Future analysis 
should focus on the relevance of corrections beyond the 
on-shell approximation, on the case A < 1, and on higher 
temperatures T > J. 

We are indebted to F. Heidrich-Meisner and R. G. 
Pereira for valuable comments. Work supported in part 
by the Deutsche Forschungsgemeinschaft through Grant 
No. BR 1084/6-1, FOR912 and by the National Science 
Foundation under Grant No. PHY05-51164. 
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